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Project  Abstract  (Original  Proposal) 

Many  experimental  investigations  have  noted  that  complex  shock  /  boundary  layer 
interactions  exhibit  large-scale,  low  frequency  motion  that  is  distinct  from  that  associated  from 
wall-layer  turbulence.  This  motion  results  in  fluctuating  pressure  and  heat-transfer  loads  that  can 
adversely  affect  the  performance  of  supersonic  missile  systems,  among  others.  Accurate 
predictions  of  such  effects  would  appear  to  require  a  time-dependent  modeling  strategy,  such  as 
large-eddy  simulation  (LES)  or  direct  numerical  simulation  (DNS).  The  application  of  these 
techniques  to  flows  at  practical  Reynolds  numbers  requires  immense  computational  resources, 
and  as  yet,  very  few  studies  have  provided  one-to-one  comparisons  with  experimental  data.  A 
promising  alternative  to  traditional  LES/DNS  are  hybrid  large-eddy  /  Reynolds-averaged  Navier- 
Stokes  (LES/RANS)  models,  which  essentially  function  as  a  large-eddy  simulation  away  from 
solid  surfaces  but  revert  to  a  RANS-based  closure  near  solid  surfaces.  Recent  simulations  of  two- 
and  three-dimensional  shock  /  boundary  layer  interactions  using  a  class  of  “zonal”  LES/RANS 
methods  developed  at  North  Carolina  State  University  have  shown  that  these  techniques  can 
predict  time-averaged  flow  field  properties  as  well  as  or  better  than  RANS  models.  The 
techniques  naturally  capture  both  large-scale  and  small-scale  flow  features,  but  an  open  question 
is  whether  the  dynamics  of  the  shock  /  boundary  layer  interactions  is  predicted  correctly.  This 
tightly  focused  research  effort  will  attempt  to  answer  this  question  by  conducting  detailed  hybrid 
LES/RANS  simulations  of  the  Mach  5,  28  degree  compression-comer  flows  studied  by  David 
Dolling’s  group  at  the  University  of  Texas.  This  data  is  the  most  comprehensive  set  available 
for  studying  shock  interaction  dynamics  and  includes  conditionally-averaged  pitot  and  static 
pressure  distributions,  wavelength,  amplitude,  and  frequency  information  associated  with  the 
shock  structure,  and  cross-correlation  analyses.  Simulations  will  be  conducted  for  this 
interaction  and  will  involve  investigation  of  the  effects  of  variations  in  the  turbulence  modeling 
parameters,  the  extent  of  the  computational  domain,  and  the  degree  of  mesh  refinement  on  the 
predictions.  Data  extraction  will  mimic  that  used  in  the  experimental  measurements  (to  the 
extent  possible).  This  will  enable  the  one-to-one,  time-resolved  comparisons  necessary  for  a 
complete  assessment.  Results  obtained  through  a  successful  validation  should  provide  directions 
for  improvement  in  less-expensive,  lower- fidelity  models  (unsteady  RANS  or  very  large  eddy 
simulation  (VLES))  and  will  perhaps  yield  new  insights  into  the  causes  of  large-scale 
unsteadiness  in  flows  driven  by  shock  /  boundary  layer  interactions. 
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1.  Background  and  Objectives 

The  accurate  prediction  of  high-speed  flows  can  depend  critically  on  the  ability  of  the 
chosen  model  to  account  for  shock  wave  /  boundary  layer  interactions.  For  example,  off-design 
aerodynamic  performance  of  missiles  can  be  affected  by  shock  /  boundary  layer  interactions 
induced  by  plume  interference,  and  the  initial  stages  of  launch  for  many  ballistics  concepts 
involve  multiple  shock  /  boundary  layer  interactions.  Maneuvering  of  launched  missiles  may 
involve  movement  of  fins  or  control  surfaces,  which  could  also  result  in  shock  /  boundary  layer 
interactions.  A  clear  understanding  of  the  effects  of  the  associated  fluctuating  pressure  and 
heating  loads  is  necessary  to  better  refine  control  techniques,  to  predict  failure  scenarios,  and  to 
determine  ranges  of  efficient  operation.  A  complicating  factor  in  many  shock  wave  /  boundary 
layer  interactions  is  a  degree  of  global  flowfield  unsteadiness,  apart  from  that  due  to  turbulent 
fluctuations,  that  can  dominate  the  response.  The  work  of  David  Dolling  and  his  colleagues  [1- 
8],  spanning  nearly  twenty  years,  has  provided  the  most  complete  characterization  of  the 
dynamics  of  shock  /  boundary  layer  interactions.  The  database  includes  conditionally-averaged 
pitot  and  static  pressure  distributions  (conditioned  on  the  separation  shock  position),  wavelength, 
amplitude,  and  frequency  information  associated  with  the  shock  structure,  and  cross-correlation 
analyses  designed  to  elucidate  the  underlying  causes  of  the  global  motion.  Trends  revealed 
during  this  study  raise  serious  concerns  about  the  ability  of  traditional  Reynolds-averaging 
techniques  to  predict  shock  /  boundary  layer  interactions  and  appear  to  indicate  that  time- 
accurate,  three-dimensional  modeling  strategies  (such  as  DNS  or  LES)  may  be  required  for 
precision. 

The  objective  of  the  current  work  is  to  apply  a  hybrid  LES/RANS  framework  developed  at 
NCSU  [9-13]  to  the  Dolling  Mach  5  compression-comer  configuration  to  determine  the  model’s 
ability  to  capture  the  time-dependent  structure  of  this  flow.  The  viewpoint  taken  in  the  NCSU 
model  is  that  the  RANS  component  functions  as  a  form  of  near-wall  closure  in  an  otherwise  LES 
description  of  the  flow.  A  strong  coupling  of  the  model  response  with  the  boundary  layer 
structure  is  therefore  required,  and  this  is  facilitated  first  by  the  design  of  a  RANS-to-LES 
transition  function  that  is  based  on  the  modeled  structure  of  the  flow  in  the  logarithmic  region 
and  second,  by  the  use  of  “recycling  /  re-scaling”  techniques  that  sustain  the  transport  of  large 
turbulent  structures  in  the  outer  part  of  the  boundary  layer.  While  predictions  of  mean-flow 
properties  have  been  evaluated  in  earlier  works  [9-13],  this  study  also  focuses  on  the  prediction 
of  second-moment  quantities  and  on  the  dynamics  of  the  shock  /  boundary  layer  interaction.  In 
this  scope  and  in  its  one-to-one  comparisons  with  available  experimental  data  at  the  same 
conditions,  the  present  work  is  closely  related  to  the  Mach  3  LES  investigation  of  flow  over  a 
compression  comer  conducted  by  Loginov,  et  al  [14],  the  Mach  3  DNS  compression-corner 
calculations  of  Wu  and  Martin  [15],  and  the  earlier  “very”  large-eddy  simulations  of  Mach  3 
flow  over  a  compression  comer  of  Hunt  and  Nixon  [16].  Both  the  Loginov,  et  al.  and  Wu  and 
Martin  investigations  used  over  18  million  mesh  points  in  their  calculations.  Because  of  the 
near- wall  modeling  in  our  hybrid  LES/RANS  approach,  we  can  reduce  the  resolution 
requirements  significantly  and  can  compute  flows  at  even  higher  Reynolds  numbers.  It  will  be 
shown  that  many  of  the  features  of  the  Dolling-group  compression  comer  experiments  are 
captured  to  a  good  degree  of  precision  by  the  hybrid  LES/RANS  model.  The  remainder  of  this 
report  outlines  the  numerical  methods  used  (Section  2),  the  turbulence  closure  models  (Section 
3),  and  other  calculation  details  (Section  4).  The  results  of  the  investigation  are  presented  in 
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Section  5,  some  conclusions  are  outlined  in  Section  6,  and  directions  for  future  work  are 
presented  in  Section  7. 

2.  Numerical  Methods 


Simulations  were  performed  using  NCSU’s  REACTMB  code  [17].  In  REACTMB,  the 
governing  equations  are  discretized  in  a  finite-volume  manner.  Inviscid  fluxes  are  discretized 
using  Edwards’  LDFSS  scheme  [18],  while  viscous  and  diffusive  fluxes  are  discretized  using 
second-order  central  differences.  The  LDFSS  scheme  is  extended  to  second  or  higher-order 
spatial  accuracy  using  the  piecewise  parabolic  method  (PPM)  [19],  as  described  later.  A  dual¬ 
time  stepping  implicit  method  is  used  to  advance  the  equations  in  time.  At  each  time  step,  a 
Crank-Nicholson  discretization  of  the  equations  is  solved  to  a  prescribed  tolerance  using  a  sub¬ 
iteration  procedure.  The  matrix  system  resulting  from  the  linearization  of  the  equation  system  is 
approximately  solved  using  a  planar  relaxation  procedure  at  each  sub-iteration.  To  enhance 
computational  efficiency,  matrix  elements  are  evaluated  and  factored  every  second  time  step  and 
are  held  fixed  over  the  duration  of  the  sub-iterations. 


In  an  effort  to  enhance  the  ability  of  REACTMB  to  capture  fine-scale  structures,  the  Sonic-A 
ENO  scheme  [20]  used  in  our  prior  work  has  been  replaced  by  a  version  of  the  classic  piecewise 
parabolic  method  (PPM)  [19].  In  contrast  to  the  ENO  schemes  described  in  [20]  and  to  the 
standard  forms  of  weighted-ENO  (WENO)  schemes,  the  PPM  reduces  to  a  fourth-order  central 
differencing  scheme  for  sufficiently  smooth  data,  enhancing  its  ability  to  capture  high-frequency 
features.  A  cell-by-cell  limiting  procedure  reduces  the  order  of  accuracy  in  the  vicinity  of  local 
extrema  in  the  solution.  A  general  interface  flux  formula  may  be  written  as 

P'i+uii^Lj+vi^R.t+i/i)  >  left-and-right  state  interpolations  VLM/2  and  VR  Mn  determined  by  the 
step-by-step  procedure  described  below. 


Step  1:  initially  set  left-and-right  state  values  to  fourth-order  averaging  operators.  If 
unmodified,  these  provide  a  fourth-order  central  difference  approximation  to  the  flux 

difference  Fi+U2  -  Fi_]l2 : 


;+l/2  —  Vr,  1+1/2  —  2 


(i) 


Step  2:  conduct  cell-by-cell  limiting  to  enforce  monotonicity: 


5 


if  sgn[  (VLMin  -KXK-  V,w)]  =  -1.  ih™ 

^i,;+l/2  =  VrJ-  1/2  =  ^ 

else 


r  =  (/  _F 

^  '  i.i+l/2  r  R,i-\!2 


D  =  6.0[V,~i(V, 
if  (DC  >  CC )  then 


Ij+1/2  +  ^,1-1/2)] 


F  =  3F  -  2F 

y  Rj-l/2  JV  i  L,  1+1/2 

elseif  (-CC  >  DC)  then 


F  =  3F 

r  L,i+U2  i 


■2V, 


R,i- 1/2 


end  if 
end  if 


(2) 


The  first  ‘if  block  resets  the  interpolation  function  to  a  constant  if  Vi  is  a  local  maximum  or 

minimum.  The  second  ‘if  block  resets  either  the  left-state  value  at  interface  i  + 1/2  or  the  right- 
state  value  at  interface  i  - 1  /  2  so  that  the  interpolation  parabola  that  connects  the  interface  states 
with  the  state  at  the  cell  center  is  monotonically  increasing  or  decreasing.  It  is  also  possible  to 
enforce  physical  constraints  (such  as  positive  temperatures,  densities,  and  mass  fractions)  at  the 
cell  interfaces  by  a  similar  cell-by-cell  resetting  algorithm.  The  PPM  requires  a  seven  point 
stencil  in  each  coordinate  direction,  and  the  reconstruction  procedures  are  applied  to  the 

primitive-variable  vector  V  =  [p,u,v,w,T,k,cof  . 


3.  Turbulence  Closure 

3.1.  Menter’s  SST  Model 


The  RANS  turbulence  model  used  in  this  investigation  is  Menter’s  hybrid  k-co/k-s  shear- 
stress  transport  (SST)  model. [21]  The  transport  equations  for  the  turbulence  kinetic  energy,  k , 
and  the  specific  dissipation  rate,®  for  Menter’s  model  are  given  by 


aW+s{phP=/i'S2_/j.paik+_d 


dt 


dx , 


aM , 


dxj 


(M  +  VkMt) 


dk 

dx. 


(3) 


dt 


-  +  - 


dx, 


dx, 


/  \  dco 

dx, 


1  dk  dco  ... 
+  2(1  -  Fx  )poal  ——  —  (4) 


co  dxj  dxj 


where  S  is  defined  as 


5  = 


dui  duj  dut  dui 


dxj  dxf  dxj  dxj 


dut  'l 

2 

W-J 

(5) 


6 


and  <f)  represents  any  constant  in  Menter’s  model  (cr k  averaged  according  to  the  blending 

function  Fl 

</>  =  Fl</>l+(l-Fl)<?>2  (6) 

The  (f>x  constants  are  from  the  k  -  co  model: 

<?k\  =  0.85,  <7ml  =  0.5  ,  J3X  =  0.075 


j3  =0.09,  /f  =  0.41,  yx  =  44~ 


A 


and  the  (j)2  constants  are  from  the  standard  k-s model: 


crk2  =1.0,  (7a2  =0.856,  P2  =0.0828 
A*  =  0.09,  k  =  o.4i ,  r2=^-^f*L 


(7) 

(8) 


The  blending  function  Fx  is  defined  as 

Fx  =  tanh(argj’ )  ,  argj  =  max 


4k  500iA 


v0.09<ai  ’  d2co  j 

where  d  is  the  distance  to  the  nearest  wall.  The  eddy  viscosity  is  defined  as: 

axk 


v,  =  ■ 


(9) 


(10) 


maxima?,  AssIQF2) 

where  a,  =  0.31 ,  Q  is  the  magnitude  of  the  vorticity  vector,  and  F2  is  another  blending  function 
given  by 

4k  500tA 


F2  =  tanh(argj),  arg2  =  max 


v  0.09 cod  d  co  j 


(ii) 


This  ‘shear  stress  transport’  (SST)  modification  alters  the  turbulence  frequency  as  used  in  the 
eddy  viscosity  definition  in  the  vicinity  of  high  strain  rates.  The  general  effect  is  to  lower  the 
eddy  viscosity  and  thus  to  promote  the  growth  of  regions  of  separated  flow.  As  shown  later,  the 
RANS  predictions  for  Dolling’s  compression-comer  experiments  (and  similar  shock  /  boundary 
layer  interactions)  are  sensitive  to  the  constant  Asst  in  Eq.  (9),  which  ranges  between  zero  (for  no 
SST  modification)  to  one  (for  the  full  SST  contribution). 


3.2.  Hybrid  LES/RANS  Extension 


Several  approaches  have  been  used  in  our  earlier  work  to  extend  Menter’s  model  to  serve  as  the 
RANS  component  of  a  hybrid  LES/RANS  closure  [9-13].  In  this  investigation,  we  have  opted 
for  a  simple  strategy  motivated  by  Strehlets’  [22]  development  of  a  two-equation  detached-eddy 
simulation  method.  The  eddy  viscosity  as  used  in  the  turbulence  transport  equations  and  the 
main  flow  equations  is  defined  as  a  weighted  sum  of  a  RANS  description  and  a  Smagorinsky 
model: 


A,  =  py,  =  p 


- - ^(i-r)CsA2/> 

max  (a2(o,  AsstklF2 )  , 


,A  =  (AAAz)1/3,Cj=0.01  (12) 


x  y 
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The  blending  function  T  is  designed  to  shift  the  closure  from  a  RANS  description  near  solid 
surfaces  to  an  LES  description  in  the  outer  parts  of  the  boundary  layer  and  in  regions  of  flow 
separation.  The  blending  function  is  based  on  the  ratio  of  the  wall  distance  d  to  a  modeled  form 
of  the  Taylor  micro-scale: 


r  = 


2 


f 


\ 


1  -  tanh[  5( 


Vc, 


d 

7  = - 

«  a 


(13) 


with  the  Taylor  micro-scale  defined  as 

X  =  ^Jc^o 


(14) 


The  primary  advantage  of  this  form  is  that  the  location  of  the  RANS/LES  juncture  (defined  as 
f  =0.5)  can  be  correlated  as  a  function  of  the  wall  coordinate  d 1  =  uTd  /  v  .  Substituting  the  log- 


law  expression  for  the  turbulence  frequency  co  =  uT  into  (14)  and  placing  the  result  in 

(13),  it  can  be  shown  that  the  r=0.5  position  should  occur  at  d+  =  a2  if  <f>=  0.0.  The  constant 
ax  then  can  be  used  control  the  position  of  the  juncture,  assuming  that  the  log-law  scaling  is 
maintained.  The  value  of  “5”  in  Eq.  (13)  controls  the  sharpness  of  the  transition,  and  the 
constant^  can  be  used  to  shift  the  mid-point  of  the  blending  function.  In  this  work,  <f>  is  set  to 


2.2975599,  so  that  the  balancing  position  (where 


VC 


V 


=  1 )  corresponds  to  T  =0.9  instead  of 


r=0.5. 


To  determine  the  constant  al  for  a  particular  inflow  boundary  layer,  the  following  procedure  is 

used.  First,  a  prediction  of  the  equilibrium  boundary  layer  is  obtained,  given  free-stream 

properties,  a  specified  wall  condition  (adiabatic  or  isothermal)  and  a  value  for  the  boundary  layer 

thickness  (0.0175  m  for  the  Dolling  Mach  5  experiments),  from  Coles’  Law  of  the  Wall  /  Wake 

along  with  the  Van  Driest  transformation: 

u.  1  ,  ,  „n  .  2,71  d.  ud 

^  =  -ln(0  +  C  +  2-sin2(--),  dl=-*~ 

UT  K  K  2  O  V 

T  W 

with 

2A2u/u„  -  B 


(15) 


Kd  = 


A  = 


A 


sin 


a/52+4^2 


+  sin 


B 


^B2  +  4A2 


(16) 


^  Pr  M2  ,  B  = 

1  +  Pr//2  (^  !)M2 

I  2  K 

2 

2^_i 


An  initial  estimate  for  the  outer  extent  of  the  log  layer  is  defined  by  finding  the  value  of  d^  such 


that 


-NO  +  C 


^  fu  \ 

vd  1  =  0.98 .  The  value  of  d+  =  u r d  !v  that  corresponds  to  this  value  of 


\Ur  J 


d+K  is  then  found  through  the  use  of  Walz’s  formula  for  the  static  temperature  distribution  within 
the  boundary  layer  [23]: 


8 


(17) 


T_ 

T 

oo 


(Taw-TJ  u  (r- 1) 


m: 


\UC0j 


Since  the  kinematic  viscosity  v  is  a  function  of  temperature,  the  target  value  for  d+  will  differ 
significantly  from  d +  for  high  Mach  number  flows.  For  the  inflow  boundary  layer  in  the  Mach 
5  Dolling  experiments,  a  value  of  ax  =42.042  is  calculated  using  this  procedure. 

A  key  to  the  blending  function  response  is  the  relative  insensitivity  of  the  turbulence  frequency 
to  changes  in  the  turbulence  kinetic  energy  and  the  eddy  viscosity.  This  means  that  the  average 
location  of  the  blending  function  (in  terms  of  its  distance  from  the  wall)  will  not  vary 
significantly.  The  dependence  of  the  turbulence  frequency  on  the  local  strain  rate  S  means  that 
the  instantaneous  position  of  the  interface  will  vary  in  response  to  resolved-eddy  dynamics.  The 
calibration  procedure  described  above  is  problem-dependent,  and  it  is  not  clear  as  yet  that  the 
procedure  is  self-similar  in  that  the  T  =0.5  blending  function  position  would  recover  to  its  target 
location  if  the  initial  equilibrium  boundary  layer  relaxes  to  another  equilibrium  state.  Efforts  are 
underway  to  analyze  the  response  of  this  technique  (and  other  blending  strategies)  using  RANS 
mean  data. 


3.3.  Recycling-Re-scaling  Method 


A  recycling  /  rescaling  technique,  applied  to  the  fluctuating  fields,  is  used  to  initiate  and  sustain 
turbulent  structures.  In  this  procedure,  fluctuations  in  the  fluid  properties  at  a  ‘recycle  plane’  are 
extracted  by  subtracting  the  instantaneous  profile  from  a  time-  and  span-averaged  profile.  The 
fluctuation  fields  are  then  rescaled  according  to  boundary  layer  similarity  laws  and  superimposed 
onto  a  RANS  mean  inflow  profile.  The  general  procedure  is  described  in  [10]  and  [11].  One 
modification  used  in  this  work  to  prevent  excessive  turbulence  energy  accumulation  in  the  outer 
part  of  the  boundary  layer  is  to  multiply  the  recycled  fluctuations  by  a  Klebanoff-type 
intermittency  function  [24]: 

inflow  =  7rans  +  ^kieb  7,'ecyc  J  7  =  density,  velocity,  temperature 


inflow  =  (1  -  ^kieb  )7rans  +  ^kieb  W  J  7  =  turbulence  variables  (18) 


Faeb=(i+(d/cklebs)6y 

with  Ckleb=  1.10.  This  procedure  also  ensures  that  the  RANS  free-stream  inflow  properties  are 

not  altered  outside  the  boundary  layer.  In  our  prior  work,  we  have  assumed  negligible  pressure 
fluctuations  to  determine  the  fluctuations  in  density,  given  the  fluctuations  in  temperature.  We 
relax  this  assumption  in  this  work  in  order  to  allow  a  user-specified  pressure  fluctuation  level. 
We  also  require  that  the  recycled  temperature  fluctuations  be  no  greater  (in  magnitude)  than 
allowed  by  Morkovin’s  hypothesis  of  negligible  total-temperature  fluctuations.  From 
Morkovin’s  hypothesis, 


C  T'  = 

p 1  mork 


kleb  (^RANS  Wrecyc  +  VRANSVrecyc  +  ^RANS^ 


recyc 


)  + 


2  ^kleb  (^ recyc  Fecyc  Aecyc  ) ) 


(19) 


and  the  recycled  temperature  fluctuation  is  limited  as  follows: 
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Figure  1 :  Isosurfaces  of  vorticity  magnitude  colored  by  axial  velocity 


T' 

recHlim 


=  ^klebCyc  mm(]  A 


V 

mork 


F  T' 

kleb  recyc 


), 


^kleb^recyc  ^  0 


(20) 


The  density  fluctuation  is  determined  as  follows.  First,  a  provisional  value  of  the  pressure 
fluctuation  is  determined  from  the  recycled  density  fluctuation: 


prov 


^ I- /^RAN S ^recyc  |  j ^kleb Pvqcjc  ( ^RAN S  ^recyc|l;,„  )1 


rec^c  I  lim  J 


(21) 


The  pressure  fluctuation  is  then  limited  to  be  a  specified  multiple  of  the  pressure  (2%  in  this 
case): 


^ProvL  =Sgn(/?prov)min(|^prov|5°-02^RANs)  (22) 


and  a  corrected  value  of  the  density  fluctuation  is  determined  from  the  limited  pressure 
fluctuation: 


reeve  =  (P  prov  I  lim  [  )  /(Trans  +  ^recyc  |  ,:m  1 


lim 


recyc  I  lim  } 


(23) 


4.  Calculation  Details 

4.1.  Test  Geometry  and  Boundary  Conditions 

The  test  geometry  is  a  28  degree  compression  ramp  placed  on  the  floor  of  the  University  of 
Texas  Mach  5  blowdown  wind  tunnel. [5]  Several  specific  ramp  configurations,  differing 
primarily  in  the  width  of  the  ramp,  have  been  used.  The  one  modeled  in  this  investigation 
(Ramp  3  in  [5])  is  3.5  inches  wide.  Aerodynamic  fences  are  attached  to  both  sides  of  the  ramp. 
The  precise  dimensions  of  the  fences  are  not  given  in  [5],  though  it  is  stated  that  the  leading  edge 
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of  the  fence  at  the  tunnel  wall  extends  2.75  inches  upstream  of  the  ramp  leading  edge.  Figure  1 
shows  isosurfaces  of  vorticity  magnitude  (50000  s'1),  colored  by  axial  velocity.  The  position  of 
one  of  the  fences  is  indicated.  The  effects  of  the  fences  are  included  in  the  three-dimensional 
calculations  by  imposing  solid  wall  boundary  conditions  over  the  portions  of  the  Zmm  and  Zmax 
grid  planes  that  are  within  the  fence  area.  Otherwise,  periodic  boundary  conditions  are  imposed 
at  these  planes.  The  mesh  spacing  in  the  spanwise  (Z)  direction  is  uniform,  thus  boundary  layer 
growth  on  the  fences  is  not  resolved.  The  tunnel  /  ramp  surface  (Ymin)  is  treated  as  an  adiabatic 
wall.  An  extrapolation  condition  is  used  at  the  upper  boundary  (Ymax)  and  at  the  outflow 
boundary  (Xmax).  The  recycling  /rescaling  technique  mentioned  above  is  used  to  provide  a  time- 
dependent  inflow  boundary  condition  at  Xmin.  The  computational  domain  extends  from  Xmm  = 
0  m  to  Xmax  =  0.304  m  in  the  streamwise  direction,  from  Zmm  =  -0.0444  m  to  Zmax  =  0.0444  in  the 
spanwise  direction,  and  from  Ymjn  =  0  m  to  Ymax  =  0.0635  m  (relative  to  the  plate  surface)  in  the 
vertical  direction.  The  ramp  leading  edge  is  located  at  X  =  0.1905  m.  The  baseline  mesh 
contains  360  cells  in  the  X  direction,  144  cells  in  the  Y  direction,  and  129  cells  in  the  Z  direction 
(6.687  million  cells  in  total).  The  mesh  is  clustered  to  the  plate  surface  using  a  hyperbolic 
tangent  stretching  function  so  that  the  minimum  d+w  in  the  inflow  boundary  layer  is  around  one. 

The  mesh  spacing  is  uniform  in  the  X  and  Z  directions  and  is  chosen  such  that  approximately  20 
cells  would  be  located  within  one  boundary  layer  thickness  (0.0175  m).  One  simulation  was 
performed  on  a  mesh  with  smaller  spacing  (60%  of  the  baseline  value)  in  the  interaction  region. 
This  mesh  contains  480  cells  in  the  X  direction,  thus  raising  the  total  number  of  cells  to  8.916 
million.  For  the  remainder  of  this  paper,  the  X  coordinate  is  redefined  to  represent  the  distance 
along  the  surface  of  the  geometry,  and  the  origin  is  located  at  the  wedge  apex.  The  Y  coordinate 
is  redefined  to  represent  the  coordinate  normal  to  the  surface. 

4.2.  Initialization  Procedure  and  Time  Evolution  Details 

The  calculations  were  initialized  by  first  computing  a  2-D  flat  plate  solution  for  the  incoming 
boundary  layer.  A  target  inflow  profile  was  determined  as  the  location  where  the  computed 
boundary  layer  properties  best  matched  the  experimental  values  shown  in  Table  1. 

Table  1:  Inflow  boundary  layer  properties 


Parameter 

Value 

M,x 

4.95 

Po  (Pa) 

2.17e6 

T0(K) 

350 

Re/m 

50.1e6 

50  (cm) 

1.75 

8*  (cm) 

0.63 

0  (cm) 

0.00513 

cf 

7.83e-4 

The  test  section  in  the  experiments  is  just  downstream  of  the  wind-tunnel  nozzle  exit.  Details  of 
the  nozzle  geometry  were  not  available,  so  the  above  procedure  was  used  to  provide  a 
representative  inflow.  The  inflow  profile  was  extracted  from  the  flat-plate  data  and  used  to 
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initialize  a  2-D  RANS  calculation  of  the  entire  interaction.  As  discussed  in  Section  5.5,  the 
choice  of  the  constant  Asst  has  a  pronounced  effect  on  the  prediction  of  the  upstream  extent  of 

axial  separation;  the  2-D  RANS  data  corresponding  to  Asst=  0.9  was  used  to  initialize  the  three 
dimensional  hybrid  LES  /  RANS  calculations. 


To  accelerate  the  development  of  coherent  structures,  a  velocity  fluctuation  field  obtained  from  a 
previous  LES/RANS  simulation  of  a  Mach  3  boundary  layer  at  a  Reynolds  number  /  meter  of 
49.5e6  [25]  was  superimposed  upon  the  RANS  base  state  upstream  of  the  wedge  apex.  The 
Mach  3  velocity  fluctuation  field  was  expanded  in  the  periodic  Z  (spanwise)  and  assumed 
periodic  X  (streamwise)  directions  to  cover  the  X-Z  extent  of  the  flat-plate  region.  The  velocity 
fluctuations  from  the  Mach  3  simulation  were  scaled  by  the  ratio  of  the  friction  velocities, 

(24) 


r 

I  A 

"a 

ll 

uT 

ir, 

II 

\M=5  1  IM=3 

[ut 

CO 

II 

and  all  coordinates  were  scaled  by  the  ratio  of  the  boundary  layer  thicknesses, 


x, 


1  -  1 

\M= 5  Xi  \M= 3 

\ 

M= 5 

u 

M= 3  J 

(25) 


to  provide  a  reasonable  match  with  the  Mach  5  flow.  Equation  (19)  was  then  used  to  provide 
the  temperature  fluctuation  field,  and  an  assumption  of  negligible  pressure  fluctuations  was  used 
to  provide  the  density  fluctuation  field.  Wu,  et  al.  [26]  have  reported  the  use  of  a  similar 
initialization  technique  in  their  DNS  calculations. 


Starting  from  this  initial  condition,  a  hybrid  LES/RANS  solution  with  Assl  =  0.9  was  evolved  at  a 
time  step  of  2e-7  s  for  3.5  ms  (155  S0/ux)  to  eliminate  initial  transients.  Complete  time 
histories  for  several  planes  of  data  (Z J S0  =  0.0,  Y/ S0  =  0.0,  Y/ S0  =  0.2,  Y/  S0  =  0.7,  X/ S0  =-3.5, 
X/ S0  =  0.0,  X/  S0  =  4.0)  were  then  collected  for  an  additional  8  ms  (354  S0  /  ux ).  The  sampling 
period  was  2e-6  s.  Flow-property  data  was  also  averaged  in  time  and  in  space  over  the  interval 
-1.8<  z  /s0  <  1.8  to  yield  nominally  two-dimensional  statistics.  The  solution  obtained  after 

11.5  ms  was  used  as  the  initial  condition  for  three  additional  runs:  Asst  =  0.0,  Asst  =  0.9  on  the 
refined  mesh,  and  Asst=  1.0.  In  these  cases,  the  solutions  were  evolved  for  3.0  ms  without 
gathering  statistics.  Statistics  were  then  collected  for  a  minimum  of  8  ms  beyond  this  point. 

5.  Results  and  Discussion 


5.1.  Flat  Plate  Boundary  Layer 

Axial  velocity  profiles  in  wall  coordinates  are  compared  with  experimental  data  [5]  in  Figure  2. 
These  profiles,  as  well  as  those  shown  in  Figure  3,  were  extracted  at  the  recycle-plane  location 
(-7 1)0  downstream  of  the  inflow  plane).  The  RANS  and  hybrid  LES/RANS  solutions  slightly 

over-predict  the  experimental  data  but  are  virtually  indistinguishable  from  one  another.  The 
time-averaged  position  of  the  blending  function  is  also  indicated.  The  r=0.9  position  is  at 
<7^=266.1,  while  the  r=0.5  position  is  at  <7*  =344.5.  The  target  d*  location,  defined  according 
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Figure  2:  Inflow  boundary  layer  velocity  profile  and  blending  function  position  (wall 
coordinates) 

to  the  criteria  outlined  in  Section  3.2,  is  345.5,  so  it  is  clear  that  the  predicted  average  position  of 
the  blending  function  is  obtained  in  the  simulations. 

Figure  3  compares  Reynolds  stress  components  and  rms  mass  flux  fluctuation  data  with 
experimental  data  from  Smits  and  Muck  [27]  at  Mach  2.79,  Re/m  =63.0e6,  as  comparable 
Reynolds  stress  data  for  the  actual  Mach  5  experiments  could  not  be  located.  The  Reynolds 
stress  components  are  scaled  by  pjp  while  the  rms  mass  flux  fluctuation  is  scaled  by  pjix ,  thus 
providing  a  measure  of  the  turbulence  intensity.  Agreement  with  experimental  data  for  the 
normal  stress  component  <  pu  u  >  is  good,  but  the  predictions  for  the  resolved  shear  stress 
component  <  puv  >  show  a  larger  constant-stress  region  (extending  to  approximately  Y/<5loc  = 

0.4)  that  is  indicated  in  the  experimental  data.  This  leads  to  an  over-prediction  of  the  Reynolds 
shear  stress  in  the  outer  part  of  the  boundary  layer.  The  sum  of  the  modeled 
<  pu  v  >  component  (also  shown)  and  the  resolved  <  pu  v  >  component  extends  the  constant- 
stress  region  further  toward  the  wall,  and  the  maximum  normalized  shear  stress  is  around  one,  as 
expected  from  inner-layer  scaling  theory.  The  prediction  for  the  normal  stress 
component  <  pvv  >  also  displays  a  maximum  value  at  Y/  <5|0C  ~0.4,  which  is  outside  the 

logarithmic  layer.  This  unexpected  behavior  seems  to  be  related  to  the  response  of  the  recycling 
/  rescaling  method  as  applied  to  wall-normal  velocity  component  and  is  the  subject  of  current 
study.  The  mass  flux  fluctuation  profile  agrees  reasonably  well  with  the  experimental  data 
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Figure  3:  Reynolds  stress  and  mass-flux  fluctuation  profiles  in  inflow  boundary  layer 
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Figure  4:  Isosurfaces  of  axial  vorticity  (side  view) 

though  the  peak  value  is  slightly  over-predicted.  Except  for  the  Reynolds  shear  stress  response 
in  the  outer  layer  (which  could  lead  to  faster  boundary  layer  growth  that  anticipated),  the 
equilibrium  boundary  layer  structure  appears  to  be  well-predicted  by  the  hybrid  LES/RANS 
model. 

5.2.  Three-Dimensional  Flow  Structure 

Some  aspects  of  the  instantaneous  flow  structure  are  illustrated  in  Figures  1,  4  and  5.  Vorticity 
iso-surfaces  in  Figures  1  and  4  reveal  an  increase  in  the  length  scale  and  an  increase  in  the 
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coherence  of  the  turbulent  structures  as  they  move  through  the  shock  system.  Smaller-scale 
structures  seem  to  be  attenuated  as  the  flat  plate  boundary  layer  shifts  to  a  wake-like  pattern  upon 
detaching  from  the  surface,  and  upon  re-attaching  on  the  ramp,  larger  vortical  structures  are 
formed.  Figure  5  shows  an  isosurface  of  zero  axial  velocity  colored  by  the  mean  temperature. 
The  regions  of  separated  flow  are  not  continuous  across  the  span  of  the  domain  and  they  taper  to 
a  point  when  re-attaching  on  the  ramp  surface.  High  wall  temperatures  are  found  in  the  vicinity 
of  the  reattachment  locations. 

Time-averaged  predictions  of  wall  skin  friction  are  shown  in  Figure  6,  and  corresponding  line 
plots  are  shown  in  Figure  7.  Results  from  three  different  simulations,  Asst  =  0.0,  Asst  =  0.9,  and 

Asst=  1.0,  are  illustrated  in  Figure  6.  The  incoming  boundary  layer  and  the  region  of  reversed 

flow  (darker  blue)  are  nominally  two  dimensional.  Some  small  peaks  and  valleys  in  the  skin 
friction  distributions  are  evident  upstream  of  the  wedge  apex,  but  the  positions  are  different 
among  the  simulations.  This  indicates  that,  with  sufficient  time  evolution,  a  statistically  two- 
dimensional  flow  would  probably  be  obtained  upstream  of  the  wedge  apex.  The  situation  is 
different  in  the  re-attachment  /  recovery  region.  Here,  the  peaks  and  valleys  in  the  distributions 
are  much  more  prominent,  with  the  regions  of  lower  skin  friction  corresponding  to  the  extension 
of  pockets  of  low-momentum  fluid  further  onto  the  ramp  surface.  In  their  LES  simulation  of 
Mach  3  flow  over  a  compression  comer,  Loginov,  et  al.  [14]  also  observe  this  topology  and 
attribute  it  to  the  presence  of  pairs  of  counter-rotating  streamwise  vortices  in  the  re-attaching 
shear  layer.  Experimental  oil-flow  images  presented  in  [28]  also  show  evidence  of  these 
structures.  The  regions  of  lower  skin  friction  (green  in  the  figure)  correspond  to  the  action  of  the 
vortices  in  lifting  fluid  away  from  the  wall,  while  the  regions  of  higher  skin  friction  (red) 
correspond  to  the  vortices’  forcing  higher-momentum  fluid  toward  the  wall.  In  the  present 


15 


Figure  6:  Contours  of  time-averaged  wall  skin  friction 

simulations,  the  spanwise  positioning  of  the  structures  correlates  with  weaker  spanwise 
inhomogeneities  in  the  incoming  boundary  layer;  the  same  effect  is  seen  in  the  simulations  of 
Loginov,  et  al.  [14]  The  distance  between  lines  of  convergence  (green  contours)  provides  a 
measure  of  the  size  of  the  dominant  vortex  pairs  We  find  that  these  sizes  range  from  ~  1.2  S0  to 

2.0  £>0 ,  which  may  be  compared  with  a  value  of  ~2.0  S0  in  Loginov,  et  al.’s  calculations. 

Loginov,  et  al.  also  argue  that  the  growth  of  these  structures  may  result  from  a  Gortler  instability 
mechanism,  given  the  degree  of  streamline  curvature  present  as  the  boundary  layer  separates, 
then  re-attaches  on  the  wedge  surface.  Navarro-Martinez  and  Tutty  [29]  have  observed  the 
formation  and  evolution  of  Gortler-type  vortices  in  simulations  of  hypersonic  laminar  flow  over 
a  compression  ramp.  The  size  of  the  vortex  pairs  in  their  work  is  approximately  the  same  as  the 
incoming  boundary  layer  thickness.  As  the  vortical  structures  appear  to  originate  from 
amplification  of  inhomogeneities  present  in  the  inflow  boundary  layer,  the  time-averaged 
positions  of  the  vortices  on  the  ramp  surface  are  likely  to  depend  on  the  methods  used  to  initiate 
and  sustain  the  boundary  layer  and  possibly  on  the  total  time  elapsed  during  the  simulation. 


Figure  7  presents  line  plots  of  time-averaged  skin  friction  at  different  spanwise  positions  for  the 
Asst  =  0.0  LES/RANS  case.  Also  shown  is  the  span-averaged  value  and  the  result  from  a  RANS 

calculation  with  Asst=  0.87.  As  shown  later,  this  choice  for  Asst  provides  the  best  match  with  the 

upstream  position  of  the  separation  shock  as  predicted  by  the  LES/RANS  method.  The  degree 
of  spanwise  variation  about  the  mean  value  is  noteworthy  (-18%  to  +26%  at  the  peak  location), 
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Figure  7:  Spanwise  variation  of  time-averaged  skin  friction 

and  at  all  stations,  the  skin  friction  values  in  the  recovery  region  exceed  those  predicted  by  the 
RANS  model.  The  high-frequency  noise  in  the  skin- friction  distributions  on  the  wedge  surface 
is  an  artifact  of  a  numerical  precision  problem  in  the  data-reduction  code  and  is  not  a  feature  of 
the  computed  flow  field.  Other  comparisons  with  the  RANS  models  are  described  in  Section 
5.5. 


The  presence  of  a  range  of  frequencies  in  the  unsteady  motion  of  shock  wave  as  it  interacts  with 
a  boundary  layer  is  well-documented  in  experiments.  ([1]  and  references  cited  within). 
Turbulent  eddies  directly  interact  with  the  shock  wave,  producing  ripples  in  the  shock  front  that 
change  their  amplitude  and  positions  over  time  scales  representative  of  the  motion  of  the 
individual  eddies.  Less  well  understood  is  the  presence  of  a  large-scale,  low  frequency 
oscillation  of  the  entire  shock  system  and  the  associated  region  of  separated  flow.  [1]  Recent 
wide-field  planar  light-scattering  (PLS)  measurements  [8]  suggest  that  this  behavior  may  be 
associated  with  the  sustained  presence  of  elongated  streaks  of  uniform-momentum  fluid  in  the 
incoming  boundary  layers.  These  streaks  have  been  shown  to  extend  upwards  of  30  S0  in  a 

Mach  2  flow.  In  regions  where  the  streaks  have  lower  momentum  than  the  average,  the 
separation  shock  can  move  further  upstream;  the  converse  is  true  in  regions  where  the  streaks 
have  higher  momentum  than  the  average,  ft  is  argued  that  the  persistence  of  these  streaks  allows 
a  low-frequency  motion  to  develop,  while  individual-eddy  interactions  continue  to  provide  a 
high-frequency  component.  In  the  LES/RANS  simulations,  we  cannot  calculate  the  entire 
evolution  of  the  boundary  layer  in  the  wind  tunnel.  Furthermore,  the  use  of  recycling  /  rescaling 
would  appear  to  force  a  length  scale  equal  to  the  distance  between  the  inflow  plane  and  the 
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Figure  8:  Snapshots  of  axial  velocity  at  Y/8o=0.2  (0.0001  s  time  intervals) 

recycling  plane  (7  S0 )  onto  the  computed  flow.  As  shown  later,  the  frequency  associated  with 

this  length  scale  is  a  noticeable  part  of  the  spectrum  but  is  not  a  dominant  part,  and  the  technique 
can  capture  lower  frequencies  as  well  as  higher  ones.  Figure  8  presents  snapshots  of  the  axial 
velocity  within  a  plane  located  at  Y/  S0  =  0.2.  This  location  corresponds  to  that  chosen  for  wide- 

field  PLS  imaging  in  [8],  The  snapshots  are  numbered  1-10  with  each  successive  snapshot 
separated  by  0.0001  s.  Contour  levels  are  chosen  to  highlight  regions  of  lower  velocity  (dark) 
and  higher  velocity  (light)  in  the  incoming  boundary  layer.  The  separation  line  is  delineated  as 
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the  onset  of  the  dark  black  region  at  the  left  of  the  figures.  There  is  clear  evidence  of  a  streak¬ 
like  pattern  in  the  incoming  boundary  layer.  Frames  1-6  in  the  snapshot  sequence  (extracted 
from  the  Asst  =  0.9  simulation  on  the  baseline  mesh)  show  that  the  separation  region  is  located 
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Figure  10:  Snapshots  of  vertical  velocity  in  X-Y  centerplane  (0.001  s  time  intervals) 

further  upstream  toward  the  upper  side  fence  (top)  and  that  this  upstream  positioning  is 
correlated  with  a  general  darkening  of  the  contours,  indicating  the  sustained  presence  of  lower- 
momentum  fluid.  At  later  times,  higher-momentum  streaks  of  fluid  eventually  enter  this  region 
and  act  to  push  the  separation  line  further  downstream.  The  separation  line  position  near  the 
lower  side  plane  does  not  migrate  as  much  over  this  time  interval.  The  passage  of  individual 
streaks  of  fluid  into  the  separation  zone  results  in  local  perturbations  to  the  separation  line,  but  it 
is  the  collective  effect  of  several  neighboring  streaks  with  approximately  uniform  momentum 
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Figure  11:  Wall  pressure  distributions  (LES/RANS  models) 

that  appears  to  induce  the  larger-scale  dynamics.  Animations  reveal  this  process  in  more  detail 
and  are  included  as  supplements  to  this  report.  Figure  9  presents  snapshots  of  the  velocity  field 
averaged  over  successive  periods  of  0.001  s.  Frame  2  in  this  sequence  corresponds  to  the  event 
illustrated  in  detail  in  Figure  8.  Frames  1,  3,  and  4  indicate  the  effects  of  isolated  streaks  of  low  / 
high  momentum  fluid  in  locally  modulating  the  separation  line,  and  the  collective  effect  of  a 
group  of  higher-momentum  streaks  in  pushing  the  separation  line  further  downstream  is  shown 
at  the  bottom  of  Frame  5. 

Figure  10  presents  snapshots  of  vertical  (v)  velocity  contours  in  the  interaction  region,  extracted 
every  0.001s  along  the  Z /<50=  0.0  centerplane  from  the  Asst=  0.0  simulation  data.  The  view 

corresponds  approximately  to  that  imaged  by  Noel  Clemens’  group  at  the  University  of  Texas 
[30]  using  time-resolved  PLS  techniques.  The  snapshots  show  a  highly  dynamic  motion  of  the 
separation  shock.  Its  progression  forward  (Frames  4,  6,  and  8)  appears  to  be  accompanied  by  the 
thickening  of  the  separation  region  (indicated  in  dark  contours)  while  its  travel  upstream  (Frames 
1  and  3)  is  associated  with  the  elongation  and  thinning  of  the  separation  region.  The  boundary 
layer  changes  dramatically  as  the  fluid  moves  through  the  separation  shock  and  re-attaches  on 
the  ramp  surface.  Larger,  more  two-dimensional  eddies  form  in  the  wake-like  region  above  the 
separation  zone,  and  their  passage  downstream  deforms  the  re-attachment  shock  plane 
significantly.  The  animations  included  in  this  report  provide  more  details  of  this  process. 

5.3.  Comparisons  with  Experimental  Mean  Flow  Data 

Mean  flowfield  data  obtained  in  the  Dolling-group  experiments  includes  average  wall  pressure 
distributions  (some  conditioned  on  the  separation  shock  position),  pressure  fluctuation  variance 
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Figure  12:  Time-and  span-averaged  axial  velocity  contours  in  the  interaction  region 


distributions,  and  pitot  pressure  distributions  in  the  re-attachment  region.  Most  of  these 
measurements  were  taken  along  the  centerline  of  the  configuration.  To  generate  smoother 
statistics  and  thereby  approximate  a  longer  time  sequence  than  is  actually  available  in  the 
simulations,  we  have  elected  to  average  most  of  data  over  the  interval  — 1.8<  Z!S„  <  1.8  in  the 

spanwise  direction.  This  eliminates  regions  of  the  flow  that  are  nearest  to  the  fences  from 
directly  affecting  the  time  averages.  Wall  pressure  distributions  are  shown  in  Figure  1 1  for  each 
of  the  LES/RANS  simulations.  Several  trends  are  apparent.  While  the  general  shape  of  the 
pressure  distribution  upstream  of  the  wedge  is  predicted  accurately,  all  of  the  calculations  under¬ 
predict  the  pressure  level  in  this  region.  Since  the  separation-shock  pressure  rise  is  a  result  of  the 
displacement  effect  of  the  separation  bubble,  this  trend  indicates  that  no  approach  correctly 
predicts  the  vertical  extent  of  the  separation  region.  The  separation-shock  pressure  rise  appears 
to  take  place  gradually  over  a  distance  of  about  1  to  1.5  S0 .  This  response  is  a  result  of  the 
averaging  effects  of  the  separation  shock  motion,  and  all  predictions  capture  this  feature  to  good 
accord.  The  Asst  =  0.0  calculation  provides  the  best  overall  agreement  with  the  experimental  data 
in  the  separation  region,  even  though  the  upstream  extent  of  the  separation  bubble  is  slightly 
under-predicted.  The  Asst=  1.0  calculation  results  in  an  elongated  separation  bubble  that  provides 

only  a  small  displacement  effect.  The  Asst  =0.9  calculations  performed  on  the  baseline  and 

refined  meshes  provide  generally  similar  results,  but  it  is  clear  that  localized  mesh  refinement 
does  not  improve  agreement  with  experimental  data.  The  pressure  level  on  the  wedge  surface  is 
in  excess  of  that  measured  in  the  experiments.  This  is  a  consequence  of  the  under-prediction  of 
the  size  of  the  separation  region  and  the  separation-shock  strength.  Other  differences  between 
the  model  predictions  are  noted  in  Figure  12,  a  plot  of  time-averaged  axial  velocity  contours  in 
the  interaction  region.  The  vertical  axis  is  the  normal  distance  from  the  surface.  As  such,  there 
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Figure  13:  Time-  and  span-averaged  kinematic  eddy  viscosity  contours  in  the  interaction 
region 

is  a  discontinuity  at  the  wedge  apex,  and  all  data  is  interpolated  onto  normal  lines  extending 
away  from  the  surface.  Also  shown  in  the  comparison  is  a  prediction  from  the  baseline  RANS 
model  with  Asst=  0.87.  The  enhanced  vertical  displacement  of  the  separation  bubble  in  the 

Ass,=  0.0  LES/RANS  calculation,  relative  to  that  of  the  Asst= 1.0  calculation,  is  clearly  indicated. 
The  bubble  size  and  shape  in  the  Assr  =0.0  LES/RANS  calculation  is  comparable  to  that  predicted 

by  the  RANS  model,  but  the  velocity  field  within  the  bubble  itself  is  clearly  different. 
Kinematic  eddy  viscosity  profiles  (normalized  with  respect  to  the  ffee-stream  kinematic 
viscosity)  shown  in  Figure  13  for  the  Asst=  0.0  and  Asst  =1.0  LES/RANS  models  indicate  the 

expected  result  of  a  lower  eddy  viscosity  when  the  SST  modification  is  fully  active.  This 
appears  to  lower  the  resistance  of  the  near-wall  flow  to  the  adverse  pressure  gradient,  allowing 
more  upstream  propagation  to  take  place  but  in  turn,  reducing  the  thickness  of  the  separation 
bubble.  It  is  clear  from  the  contours  that  this  reduction  in  bubble  displacement  weakens  the 
separation  shock. 

Predictions  of  the  standard  deviation  of  the  wall  pressure  fluctuations  are  presented  in  Figure  14. 
The  general  shape  of  the  experimental  distribution  is  captured,  but  the  peak  in  fluctuation 
intensity  near  the  separation  shock  is  under-predicted  by  all  models.  The  Asst  =0.9  calculation  on 
the  baseline  mesh  best  predicts  the  position  of  the  peak.  It  should  be  noted  that  the  magnitude  of 
the  peak  (but  not  its  location)  is  well-predicted  by  the  Asst=  0.0  calculation  if  the  data  is  sampled 

at  Y/  £>0  =0.2  instead  of  at  the  wall.  This  result  may  imply  that  the  enhanced  eddy  viscosity 

provided  by  RANS  model  in  the  near- wall  region  damps  the  separation-shock  motion  to  some 
degree.  All  calculations  capture  the  fluctuation  intensity  level  in  the  pressure  plateau  region 
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Figure  14:  rms  pressure  fluctuation  distributions  (hybrid  LES/RANS  models) 


o  data 


Figure  15:  rms  pressure  fluctuation  distributions  normalized  by  local  wall  pressure 
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Figure  16:  Time-  and  span-averaged  blending  function  in  the  interaction  region  (vertical 
coordinate  is  normal  to  the  surface) 


Figure  17:  Pitot  pressure  profiles  in  re-attachment  region 


upstream  of  the  comer,  but  all  over-predict  the  intensity  level  in  the  re-attachment  and  recovery 
regions.  As  the  wall  pressure  itself  is  also  over-predicted  in  this  region,  it  is  useful  to  examine 
the  rms  values  normalized  by  the  local  wall  pressure.  Figure  15  shows  that  none  of  the 
calculations  correctly  capture  the  relative  distributions  of  fluctuation  intensity  due  to  separation- 
shock  motion  and  due  to  shear  layer  re-attachment  and  re-attachment  shock  motion.  Whereas 
the  under-prediction  of  the  peak  in  fluctuation  intensity  near  the  separation  shock  may  be  due  to 
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Figure  18:  Separation-shock  intermittency  distributions 


the  damping  effect  of  the  near- wall  RANS  model,  the  opposite  may  be  true  on  the  ramp  surface. 
Figure  16  plots  the  time-  and  span-averaged  blending  function  in  the  interaction  region. 
Downstream  of  re-attachment,  the  blending  function  collapses  toward  the  wall,  meaning  that 
more  of  the  recovering  boundary  layer  is  modeled  as  a  large-eddy  simulation.  The  near-wall 
resolution  may  be  insufficient  to  capture  the  proper  decay  of  turbulence  energy,  and  the 
fluctuation  intensity  level  at  the  surface  may  be  over-predicted  as  a  result.  The  collapse  in  the 
blending  function  is  a  direct  result  of  an  increase  in  the  turbulence  frequency  co  as  the  flow 
passes  through  the  shock  system.  Efforts  are  underway  to  understand  this  behavior  better  and  to 
develop  modifications  that  may  control  the  response.  Even  with  this  effect,  it  is  likely  that 
further  improvements  in  the  prediction  of  the  separation-shock  strength  and  the  upstream  extent 
of  the  separation  zone  would  result  in  better  predictions  of  the  fluctuation  intensity  throughout 
the  interaction. 

Predictions  of  pitot  pressure  at  several  stations  just  downstream  of  the  comer  are  compared  with 
experimental  data  from  [4]  in  Figure  17.  The  best  result  is  provided  by  the  Asst  =  0.0  calculation, 

which  captures  the  wake-like  structure  of  the  re-attaching  boundary  layer  correctly  at  all 
positions.  For  the  other  models,  the  peaks  in  the  pitot  pressure  distribution  are  positioned  closer 
to  the  wall  and  the  boundary  layer  appears  to  accelerate  more  rapidly  than  is  indicated  in  the 
experimental  data. 
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Figure  19:  Separation-shock  intermittency  and  probability  of  reversed  flow 


5.4.  Dynamics  of  Shock-System  Motion 

References  [5]  and  [6]  describe  the  use  of  conditional  sampling  techniques  to  examine  the 
motion  of  the  shock-induced  separation  region  and  the  resulting  effects  on  the  wall  pressure 
profiles.  We  have  performed  similar  analyses  using  the  wall  pressure  data  extracted  in  the 
LES/RANS  simulations.  Again,  spatial  homogeneity  is  assumed  over  the  -1.8<  Z 1 80<  1.8 

interval  in  order  to  increase  the  sample  size.  Figure  18  shows  the  intermittency  of  the 
separation-shock  position.  At  a  particular  X  station,  the  intermittency  (y)  is  defined  as  the 
fraction  of  time  that  the  separation  shock  is  ahead  of  this  position.  A  value  slightly  greater  than 
zero  corresponds  to  the  most  upstream  positioning  of  the  separation  shock,  while  a  value  slightly 
less  than  one  corresponds  to  the  most  downstream  positioning  of  the  separation  shock.  The  most 
probable  position  of  the  separation  shock  is  indicated  by  an  intermittency  value  of  0.5.  The 
general  shape  of  the  experimental  intermittency  distribution  is  captured  by  all  LES/RANS 
models,  though  the  extent  of  the  separation-shock  motion  appears  to  be  slightly  under-predicted. 
The  most  probable  position  of  the  separation  shock  in  the  experiment  is  around  -1.8  X/  S0 ; 

the^JS<=  1.0  prediction  (-1.95  X/  S0 )  and  the  Asst=  0.9  prediction  on  the  baseline  mesh  (-1.7 
X/  S0 )  bracket  the  experimental  result.  Figure  19  plots  the  separation-shock  intermittency 

distribution  along  with  a  measure  of  the  probability  of  having  reversed  axial  flow  near  the 
surface.  A  value  of  one  at  a  given  location  means  that  for  the  entire  sample  time,  the  axial 
velocity  is  negative.  As  shown,  there  is  no  station  that  maintains  reversed  flow  at  all  times,  and 
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at  all  stations  within  the  separation  zone,  there  is  at  least  a  twenty-percent  probability  of  having 
positive  axial  flow  near  the  surface.  In  the  region  of  separation-shock  intermittency,  the 
probability  of  reversed  flow  ranges  from  near  zero  at  y  =  0  to  0.8  at  y  =  1.  From  this,  it  can  be 
concluded  that  the  most  downstream  positioning  of  the  separation  shock  is  strong  enough  to 
immediately  separate  the  boundary  layer,  while  the  separation  shock  at  its  most  upstream 
location  is  weaker  and  does  not  immediately  induce  reversed  flow.  Further  downstream,  the 
probability  of  reversed  flow  decreases  to  a  value  of  about  0.5  before  increasing  again  to  a  value 
of  0.8  near  the  wedge  apex.  A  non-zero  probability  of  having  reversed  flow  exists  on  the  wedge 
surface  past  X/  S0  =2.0.  This  is  associated  with  the  effect  of  the  counter-rotating  vortices 
discussed  in  Section  5.2  in  moving  fluid  away  from  the  surface. 

Conditional  pressure  distributions  may  be  obtained  by  ‘binning’  streamwise  wall  pressure  traces 
according  to  the  position  of  the  separation  shock.  Figure  20  compares  pressure  distributions 
corresponding  to  y  =  0  to  0.05  (the  most  upstream  positioning  of  the  shock  )  and  y  =  0.95  to  1 .0 
(the  most  downstream  positioning)  versus  experimental  data.  Good  agreement  is  generally 
achieved  for  the  most  downstream  positioning  of  the  separation  shock,  with  the  Asst  =0.0 

calculation  providing  the  best  prediction.  The  Asst  =1.0  calculation  correctly  predicts  the 

separation-shock  position  for  the  y  =  0  to  0.05  range,  but  the  pressure  level  just  upstream  of  the 
wedge  apex  is  under-predicted  and  the  pressure  level  on  the  ramp  surface  is  over-predicted. 
These  trends  again  relate  to  deficiencies  in  predicting  the  displacement  effects  of  the  separation 
bubble. 
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Figure  21:  Power  spectra  at  different  streamwise  stations 


Figure  21  presents  power  spectra  obtained  at  several  X  stations  from  pressure  data  extracted  at 
Y/c)0=0.2.  The  Asst  =0.9  database  on  the  baseline  grid  is  used  in  these  calculations;  similar  results 
are  obtained  for  the  other  databases.  As  before,  the  ensemble  of  data  in  the  interval  -1.8< 
Z/<50<  1.8  is  used.  Each  power  spectrum  is  normalized  by  the  variance  of  the  pressure  signal 
and  the  frequency,  as  in  [5].  Estimated  large-eddy  frequencies  in  the  incoming  boundary  layer 
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Figure  22:  Wall  pressure  distributions  (RANS  models) 

(based  on  S0  and  urj )  are  on  the  order  of  40-50  kHz.  The  spectrum  at  X/ S()  =-3.02  shows  that 

the  frequencies  containing  most  of  the  energy  content  lie  in  the  range  from  15  kHz  to  50  kHz. 
This  is  in  good  agreement  with  results  reported  in  [5]  (Figure  40  in  that  reference).  For  X 
stations  in  the  intermittent  region,  the  dominant  frequencies  are  much  lower:  ~0.3  to  1  kHz.  This 
is  also  in  good  agreement  with  the  results  of  Gramann  (Figure  55  in  [5])  which  show  a  dominant 
band  of  low-frequency  tones  ranging  from  -0.2  kHz  to  2  kHz.  Downstream  of  the  wedge  apex, 
the  lower-frequency  tones  diminish,  and  a  higher  frequency  spectrum  emerges.  The  dominant 
frequencies  in  the  ramp  boundary  layer  lie  in  the  interval  from  -8  to  30  kHz,  which  is  in  good 
agreement  with  Gramann’s  data  (Figure  106  in  [5]).  The  shift  to  lower  frequencies  (relative  to 
those  in  the  incoming  boundary  layer)  may  be  attributed  to  the  reduced  convection  speed  of 
eddies  downstream  of  the  shock  system,  as  the  eddy  sizes  appear  to  be  similar  in  the  calculations. 
At  all  stations,  a  selection  of  tones  with  frequencies  ranging  from  5  to  7  kHz  is  seen.  These 
appear  to  coincide  with  the  recycling  frequency,  which  is  approximately  equal  to ux  /(7S0).  This 

band  of  frequencies  is  not  dominant  at  any  location  though  its  relative  importance  does  appear  to 
increase  in  the  re-attachment  /  recovery  region.  Excepting  the  presence  of  the  recycling 
frequency,  the  predictions  are  in  close  agreement  with  experimental  observations,  indicating  that 
the  hybrid  LES/RANS  model  is  capable  of  capturing  both  low  and  high-frequency  dynamics  of 
the  interaction. 

5.5.  Comparisons  with  RANS  Mean  Flow  Data 

It  is  useful  to  compare  predictions  obtained  by  the  baseline  RANS  model  with  those  provided  by 
the  LES/RANS  model.  If  the  LES/RANS  model  can  be  proven  more  accurate,  then  it  may  be 
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Figure  23:  Skin  friction  distributions  (RANS  models) 


Figure  24:  Pitot  pressure  profiles  in  re-attachment  region  (RANS  models) 
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Figure  25:  Axial  velocity  profiles  throughout  the  interaction 


possible  to  use  LES/RANS  data  to  evaluate  modeled  forms  for  various  closure  terms  appearing 
in  RANS  models  and  from  that,  hopefully  improve  their  predictive  capability.  Figure  22  shows 
the  sensitivity  of  the  wall  pressure  distributions  predicted  by  the  RANS  model  to  the  choice  of 
the  scaling  constant^.  As  noted,  a  value  of Assl  =  0.0  results  in  the  appearance  of  a  small 

separation  bubble,  while  the  choice  of  Assl =1.0  results  in  a  very  large  separation  region.  This 
level  of  sensitivity  is  much  larger  than  that  exhibited  by  the  LES/RANS  model  (Figure  9).  The 
best  agreement  with  experimental  data  is  achieved  when  Asst  is  set  to  0.9,  while  a  value  of 

Asst=  0.87  provides  a  good  comparison  with  the  best  LES/RANS  prediction  ( Assl  =0.0).  Skin 

friction  distributions  in  Figure  23  indicate  that  the  RANS  models  consistently  predict  larger  (in 
magnitude)  near-wall  velocities  in  the  back  flow  region  in  comparison  with  the  LES  model.  The 
RANS  solutions  also  exhibit  a  delayed  re-attachment  and  a  slower  rate  of  recovery  in  the  inner 
layer,  compared  with  the  LES/RANS  solution.  Pitot  pressure  profiles  shown  in  Figure  24 
indicate  that  none  of  the  RANS  models  predicts  the  experimental  trends  as  closely  as  the 
LES/RANS  model  with  Asst=  0.0.  The  best  agreement,  not  coincidentally,  is  provided  by  the 

Asst=  0.87  RANS  solution.  Axial  velocity  profiles  at  different  streamwise  stations  in  Figure  25 
reveal  relatively  close  agreement  between  the  Asst  =0.87  RANS  solution  and  the  Asst  =0.0 
LES/RANS  solution  upstream  of  the  wedge  apex.  However,  the  transport  of  momentum  toward 
the  surface  in  the  re-attachment  region  (X/S0  =2.0)  is  clearly  larger  for  the  LES/RANS  solution. 

The  velocity  profile  further  downstream  (X/  S0  =7.0)  is  fuller  and  flatter  as  a  result,  and  the 
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Figure  26:  rms  mass  flux  fluctuation  profiles  throughout  the  interaction 

reduction  in  velocity  due  to  wall  shear  is  confined  to  a  very  thin  layer  near  the  surface.  The 
velocity  field  does  not  return  to  an  equilibrium  profile  within  the  extent  of  the  domain  for  either 
model.  In  the  absence  of  experimental  data  downstream  of  X/c)0=2.0,  it  is  difficult  to  judge 

which  approach  might  provide  the  more  correct  response,  but  other  experimental  studies  at  Mach 
3  [30,31]  have  measured  fuller,  flatter  boundary  layer  profiles  in  the  recovery  region. 
Moreover,  earlier  hybrid  LES/RANS  calculations  [10-12]  have  been  able  to  capture  the 
experimental  trends,  whereas  RANS  models  have  not. 

5.6.  Comparisons  with  RANS  Reynolds  Stress  Data 

Figure  26  presents  the  evolution  of  rms  mass  flux  fluctuation  ^[(/?«/]2)  profiles  throughout  the 

interaction  region.  Fluctuations  are  amplified  by  about  a  factor  of  10  near  the  re-attachment 
position  before  reducing  to  a  factor  of  about  7  near  the  end  of  the  ramp.  Maximum  mass  flux 
amplification  factors  of  about  4  are  reported  for  the  Mach  3  experiments  of  Smits  and  Muck  [27] 
and  Zheltovodov,  et  al.  [31]  and  for  the  computations  of  Loginov,  et  al.  [14], 

Contour  plots  of  (pk) ,  where  k  is  the  turbulence  kinetic  energy  per  unit  mass,  are  shown  in 
Figure  27  for  the  LES/RANS  model  with  Asst  =0.0  and  the  RANS  model  with  /4ss(=0.87.  Only 
the  resolved-scale  part  of  (pk)  is  presented  for  the  LES/RANS  model.  The  plots  shows  that 
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Figure  27:  Contours  of  turbulence  kinetic  energy  (normalized  by  pxii*) 

amplification  of  the  turbulence  kinetic  energy  begins  upstream  of  the  wedge,  with  the  maximum 
value  in  this  region  located  within  the  shear  layer  that  develops  above  the  separation  bubble. 
The  turbulence  amplification  effect  continues  as  the  shear  layer  re-attaches,  and  the  turbulence 
kinetic  energy  reaches  its  maximum  value  in  the  vicinity  of  the  re-attachment  shock.  Clearly 
indicated  in  the  LES/RANS  results  is  the  effect  of  the  shock  motion  in  enhancing  the 
amplification  rates  and  in  broadening  the  turbulence  kinetic  energy  distribution.  Further 
downstream,  the  turbulence  kinetic  energy  begins  to  decay,  and  the  peak  values  near  the  end  of 
the  domain  for  the  LES/RANS  and  RANS  models  are  comparable.  The  distribution  of 
turbulence  kinetic  energy  in  the  normal  direction  is  much  wider  for  the  LES/RANS  calculation, 
and  the  peak  is  located  well  away  from  the  surface.  This  may  be  due  to  convection  of 
turbulence  kinetic  energy  into  the  outer  part  of  the  recovering  boundary  layer,  as  the  mean 
velocity  gradient  in  this  region  is  relatively  small  (Figure  25).  The  RANS  result  displays  the 
expected  trend  of  locating  the  maximum  turbulence  kinetic  energy  nearer  the  wall  where  the 
predicted  mean  velocity  gradient  is  larger.  The  maximum  amplification  factor  (relative  to  the 
turbulence  kinetic  energy  level  in  the  middle  part  of  incoming  boundary  layer)  is  around  67  for 
the  LES/RANS  calculation  and  around  48  for  the  RANS  calculation.  The  amplification  factors 
reduce  to  20  and  18  for  the  LES/RANS  and  RANS  calculations,  respectively,  toward  the  end  of 
the  ramp.  These  values  are  much  larger  than  the  eight-fold  amplification  observed  in  the  Mach  3 
experiments  of  Smits  and  Muck  [27]  for  the  Reynolds  normal  stress.  In  the  Mach  3  experiments, 
the  compression  fan  that  forms  downstream  of  the  re-attachment  position  does  not  coalesce  into 
a  shock  within  the  boundary  layer,  whereas  in  the  Mach  5  calculations  (and  probably  in  the 
experiment),  the  coalescence  of  the  compression  fan  into  a  shock  takes  place  close  to  the  wall, 
and  the  coalesced  shock  interacts  dynamically  with  the  re-attaching  shear  layer. 
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Figure  28:  Contours  of  Reynolds  shear  stress  (normalized  by  ppu^) 


Contours  of  the  Reynolds  shear  stress  (pwV)  are  shown  in  Figure  28.  The  Reynolds  stress  tensor 

is  rotated  so  that  the  shear  component  is  tangential  to  the  surface.  Upstream  of  the  wedge  apex, 
the  LES/RANS  model  predicts  the  maximum  value  of  the  shear  stress  as  occurring  toward  the 
outer  edge  of  the  separating  shear  layer,  whereas  the  separation  region  is  associated  with  small 
values  of  the  Reynolds  shear  stress.  The  RANS  model  displays  the  expected  trend  of  a 
maximum  Reynolds  shear  stress  where  the  velocity  gradient  is  largest.  Downstream  of  the 
wedge  apex,  the  trends  generally  follow  that  evidenced  for  the  turbulence  kinetic  energy.  The 
maximum  amplification  factors  are  even  larger:  95  for  the  LES/RANS  calculation  and  55  for  the 
RANS  calculation.  These  reduce  to  factors  of  12.5  and  17.5  for  the  LES/RANS  and  RANS 
calculations  near  the  end  of  the  ramp.  Loginov,  et  al.  [14]  report  a  maximum  Reynolds  shear 
stress  amplification  of  28  in  their  calculations.  The  structure  factor  (pu'v'\/(pk) ,  evaluated  at  the 

location  of  maximum  Reynolds  shear  stress  toward  the  end  of  the  ramp,  is  about  0.265  for  both 
models,  indicating  that  both  computed  boundary  layers  may  be  returning  toward  equilibrium  but 
through  different  paths.  The  locations  of  the  maximum  Reynolds  shear  stress  differ 
significantly,  with  the  LES/RANS  calculation  placing  this  location  further  away  from  the  wall 
than  the  RANS  calculation.  Additional  experimental  data  would  be  required  to  determine  which 
modeling  approach  best  captures  the  evolution  of  the  Reynolds  stress  tensor  for  this  flow. 


6.  Conclusions 


Simulations  of  a  Mach  5  compression-comer  interaction  at  Re/m  =  50.1e6  have  been  performed 
using  a  hybrid  large-eddy  /  Reynolds  averaged  Navier-Stokes  (LES/RANS)  model.  The  model 
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uses  a  flow-dependent  blending  function  to  shift  the  closure  from  Menter’s  two-equation  RANS 
model  near  walls  to  a  Smagorinsky  subgrid  closure  model  in  the  outer  part  of  the  boundary  layer. 
The  blending  function  is  based  on  a  modeled  form  of  the  Taylor  microscale  and  is  designed  to 
transition  to  LES  toward  the  outer  part  of  the  logarithmic  region  (for  equilibrium  boundary 
layers).  A  recycling  /  rescaling  technique,  applied  to  the  fluctuating  fields,  is  used  to  initiate  and 
sustain  large-eddy  motion.  The  following  conclusions  may  be  stated. 

1.  The  LES/RANS  models  capture  the  mean- flow  and  time-dependent  structure  of  the 
interaction  reasonably  well.  Observed  deficiencies  relate  primarily  to  an  under¬ 
prediction  of  the  size  of  the  separation  region  (particularly  its  vertical  extent),  which 
causes  a  smaller  separation-shock  pressure  rise  and  a  larger  re-attachment  shock  pressure 
rise  than  observed  experimentally. 

2.  The  LES/RANS  results  are  sensitive  to  the  amount  of  ‘shear  stress  transport’  (SST) 
included  in  the  Menter  model  -  the  calculations  performed  without  the  SST  modification 
are  in  better  agreement  with  experimental  data.  Inclusion  of  the  SST  effect  improves 
predictions  of  the  upstream  influence  of  the  shock  system  but  reduces  the  vertical 
displacement  of  the  separation  region.  The  separation  shock  pressure  rise  is  lessened  and 
its  dynamics  not  as  vigorous  as  a  result. 

3.  In  accord  with  recent  LES  results  at  Mach  3  [14],  the  current  work  predicts  significant 
three-dimensional  effects  in  the  time-averaged  flow  downstream  of  the  re-attachment 
position.  These  appear  to  relate  to  the  formation  of  counter-rotating  vortex  pairs  in  the  re¬ 
attaching  boundary  layer.  The  positions  of  these  vortical  structures  do  not  migrate 
substantially  over  the  time  intervals  considered. 

4.  Recent  conjectures  [8]  relating  to  the  development  of  a  low-frequency  motion  of  the 
entire  shock  system  due  to  the  persistence  of  streaks  of  fluid  with  lower  or  higher 
momentum  than  the  average  are  supported  by  the  current  simulations.  In  the 
simulations,  it  appears  to  be  the  collective  effect  of  several  neighboring  streaks  that 
induces  a  low-frequency  undulation  in  the  separation-shock  front.  The  passage  of 
individual  streaks  of  fluid  through  the  front  produces  localized  ripples  but  not  large-scale 
motion. 

5.  The  dynamics  of  the  shock-system  motion,  as  quantified  by  conditional  pressure 
distributions,  shock-position  intermittency  distributions,  and  power  spectra  analysis  at 
different  locations,  have  been  predicted  to  good  accord  by  the  LES/RANS  model.  The 
use  of  a  RANS  component  near  the  surface  does  appear  to  mitigate  the  near-wall 
dynamics  of  the  separation  shock  to  some  degree.  Upon  re-attachment,  however,  the 
RANS  region  is  forced  closer  to  the  wall,  and  there  are  indications  that  the  model  may 
not  provide  enough  damping  of  large-scale  fluctuations  in  this  region. 

6.  Calculations  performed  using  the  baseline  Menter  RANS  model  are  much  more  sensitive 
to  the  amount  of  ‘shear  stress  transport’  included.  The  RANS  calculations  also  predict 
re-attachment  as  occurring  further  downstream  than  seen  in  the  experiment  or  in  the 
LES/RANS  predictions,  and  the  rate  of  recovery  of  the  boundary  layer  in  the  near- wall 
region  is  not  as  rapid  as  observed  in  the  LES/  RANS  models  or  in  the  experimental 
results  available.  Both  models  predict  substantial  amplification  of  Reynolds  normal  and 
shear  stresses  downstream  of  re-attachment,  but  the  distribution  of  the  Reynolds  stresses 
within  the  boundary  layer  is  much  broader  for  the  LES/RANS  model.  This  effect 
appears  to  be  due  to  re-attachment  shock  motion. 
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7.  From  a  perspective  of  improving  RANS  models  for  this  type  of  flow,  it  is  clear  that 
nominally  two-dimensional  flow  might  not  be  attained  downstream  of  the  re-attachment 
line.  This  means  that  two-dimensional  RANS  models  must  account  for  the  collective 
effects  of  pairs  of  counter-rotating  vortices,  and  it  is  unclear  what  average  state  one  might 
expect  to  predict.  This  fact  might  explain  why  RANS  models,  even  if  tuned  to  capture 
the  separation-shock  location,  inevitably  fail  to  capture  the  correct  recovery  rate  of  the 
boundary  layer  downstream  of  re-attachment. 

7.  Directions  for  Future  Work 

There  are  several  remaining  sources  of  uncertainty  in  the  application  of  hybrid  LES/RANS 
models  to  strongly-interacting,  wall-bounded  flows.  These  comprise  directions  for  future  work 
and  are  summarized  as  follows. 

1.  The  over-prediction  of  the  Reynolds  stress  in  the  outer  part  of  the  incoming  boundary 
layer  (Figure  3)  is  a  cause  for  concern,  assuming  that  the  Smits  /  Muck  database  is 
accurate  enough  to  be  used  as  a  benchmark.  We  believe  that  this  effect  might  result  from 
a  small  inconsistency  in  how  the  recycling  /  rescaling  method  is  being  applied,  and  we 
will  test  some  modifications  designed  to  correct  this. 

2.  The  collapse  of  the  blending  function  (and  thus  the  removal  of  much  of  the  RANS 
component  of  the  closure)  downstream  of  re-attachment  is  also  an  issue  in  that  the 
expected  dissipative  effects  of  near-wall  eddies  may  not  be  modeled  properly.  We  have 
also  observed  this  effect  in  simulations  of  a  Mach  8  crossing-shock  interaction.  [13].  We 
are  modeling  the  response  of  the  blending  functions  and  LES/RANS  turbulence  models 
using  RANS  mean  data  in  an  attempt  to  determine  how  these  models  return  to 
equilibrium.  Insights  gained  in  this  study  may  point  toward  improvements  in  the 
blending  function  approach. 

3.  The  impact  of  the  numerical  discretization  on  the  separation-shock  position  is  uncertain. 
Calculations  (not  shown)  using  our  earlier  Sonic-A  ENO  schemes  inevitably  over¬ 
predicted  the  extent  of  the  separation  region,  while  the  PPM  calculations  presented  herein 
slightly  under-predicted  the  size  of  the  separation  region.  The  capturing  of  the 
separation-shock  pressure  rise  might  be  improved  by  modifying  the  limiting  procedure  in 
Eq.  2. 

4.  We  have  not  yet  attempted  to  use  a  more  sophisticated  subgrid  closure  in  these 
calculations.  Some  candidates  include  structure-function  models  and  localized  dynamic 
models,  as  well  as  models  that  rely  on  an  equation  for  the  subgrid  turbulence  kinetic 
energy.  Our  earlier  work  used  Yoshizawa’s  model  [33]  in  this  context,  but  we  really 
haven’t  seen  a  definite  benefit  from  using  a  one-equation  subgrid  model  versus  an 
algebraic  model. 

5.  It  appears  that  predicting  the  structure  of  the  flow  in  the  re-attachment  /  recovery  region 
may  be  the  most  uncertain  aspect.  It  would  be  useful  to  re-consider  experimental 
databases  that  provide  more  data  in  this  region.  Some  candidates  include  the  Mach  3 
Smits  and  Muck  database  (and  the  earlier  work  by  Settles  and  Bogdanoff  [31]).  The 
pioneering  study  of  Dolling  and  Murphy  [34]  also  corresponds  to  these  compression- 
comer  experiments,  so  there  is  data  available  on  shock-system  dynamics. 
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6.  An  investigation  of  shock  motion  effects  as  induced  by  three-dimensional  interactions 
(sharp  fins)  may  also  be  worthwhile. 
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Appendices 

The  following  material  is  included  as  a  supplement  to  this  report. 

1. )  PDF  file  containing  John  Boles’  Master  of  Science  thesis 

2. )  Powerpoint  file  containing  a  presentation  given  at  ARO  on  May  4,  2007,  along  with 
associated  animations. 
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